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Abstract. In the recent years, several polynomial algorithms of a dy- 
namical nature have been proposed to address the graph isomorphism 
problem C |GN 02 j . |S,IC 03| . etc.). In this paper we propose a general- 
ization of an approach exposed in |GN 02| and find that this dynamical 
algorithm is covered by a combinatorial approach. It is possible to in- 
fer that polynomial dynamical algorithms addressing graph isomorphism 
are covered by suitable polynomial combinatorial approaches and thus 
are tackled by the same weaknesses as the last ones. 

1 Introduction 

In their paper, Gudkov and Nussinov proposed a new approach to analyze 
graphs. Suppose a graph has n vertices. Then each vertex of the graph is viewed 
as a point-mass in an n — 1 - dimensional Euclidean space. The point-set is 
arranged into a symmetrical initial configuration and afterward interacted by 
mutual attraction and repulsion forces, with attraction force acting on every 
pair of point-masses corresponding to vertices having an edge between them. So 
the point-set is subjected to distortion. It is possible to simulate the dynamics 
numerically, thus computing the future coordinates of the point-masses. Given 
two graphs, each of them is embedded into separate n — 1 - dimensional space. 
If after some time the two point-sets are no longer congruent, the two graphs 
are not isomorphic. However, testing the congruence of two point-sets in a d- 
dimensional space is an open problem as well, being at least as hard as graph 
isomorphism itself |A 98j . Eventually the authors have chosen the following ap- 
proach. A set of n (n — l)/2 numbers for each graph is assigned being a set of 
mutual distances between the points in the relevant point-set at a certain mo- 
ment of time. If the distance sets for the two graphs do not match, the two 
graphs are not isomorphic. The authors also conjectured that if the distances 
do match, the graphs are isomorphic. However, as pointed out in jS.TC 03| . the 
distances are the same for two non-isomorphic strongly regular graphs with the 
same parameters. 
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Essentially Gudkov and Nussinov apply their approach to give an answer to 
the graph coding problem, i.e., find a polynomial time algorithm which assigns 
an integer number to each graph, so that two graphs share the same number if 
and only if they are isomorphic. However, graph coding is not necessary to solve 
graph isomorphism, as graph vertex classification up to automorphism partition 
classes is sufficient |RC 77\ . Informally, the aim of vertex classification is to assign 
a local code to each vertex in the graph. The codes for two vertices in the same 
graph would be the same if and only if the two vertices play an equivalent role 
in the graph. 

In this paper, we shall use the idea of dynamical evolution to give a graph 
vertex classification algorithm, hereinafter referred as algorithm Al. Informally, 
given a graph on n vertices, it is embedded into n-dimensional space. Now the 
initial symmetric configuration is defined by an n x n identity matrix, which 
rows form the coordinates of the points. Afterward the point-set is subjected 
to a dynamical evolution determined by an adequate definition of attraction 
and repulsion forces acting on point-masses. If at some moment of time the 
coordinates of two points in the system are not permutation equivalent, the 
corresponding vertices are not in the same class of automorphism partition. 

In |S.IC 03j the authors show that the original |GN 02| algorithm does not 
work on strongly regular graphs and explain that by the combinatorial properties 
of strongly regular graphs. On the other hand, the algorithm Al works on many 
pairs of strongly regular graphs, including those exposed in |SJC 03j . However, 
it is still possible to indicate instances of graphs where the Al algorithm fails. 

In Section [21 we introduce the definitions and prove the basic theorems used 
in this paper. Sectional presents the dynamical system which forms the basis for 
the partitioning algorithm Al. Section 0] presents the algorithm Al, establishes 
its complexity, shows that is is covered by a polynomial combinatorial algorithm 
and indicates a counterexample, where the algorithm fails to distinguish between 
two non-isomorphic graphs. 

The fact that polynomial dynamical algorithms (at least in the presented 
case) are covered by polynomial combinatorial algorithms is the main result of 
this paper. 

In the final Section [3] it is shown that Al may be extended to a yes-no-don't- 
know type algorithm A2, which either gives a correct answer or says nothing. 

2 Preliminaries 

Definition 2.1. An (undirected) graph is an ordered pair of disjoint sets (V,r), 
where r C V x V is an anti-reflexive and symmetric relation. Elements of the 
set V are called vertices of the graph. Ifvi, v-i G V, viTv2 and v^rvi, then a set 
{«i,i>2} is called an edge of the graph. 

Thus every graph can be denoted as an ordered pair of disjoint sets G — 
(V, E), where V is the set of vertices and E - the set of edges of the graph. We 
denote by \G\ the number of vertices in the graph G. 

Graphs G\ — (Vi,Ei) and G 2 = (V2, E 2 ) are called separate, if V\ D V2 = 0- 



Definition 2.2. An adjacency matrix of a graph G = (V, E) on n vertices is 
an n x n 0,1-matrix H, such that hi t j = 1 iff {vi,Vj} G E. 

Definition 2.3. A path on a graph G — {V,E) is a sequence v\, Va, • • • , v m , 
Vi G V, such that {vx, V2}, {^2, ^3}, . . . , {v m -i, v m } are edges of the graph. 

Definition 2.4. Two vertices t>i,t>2 are called l-connected in a graph G, if ex- 
ists a path Vi,...,V2 of I + 1 elements on the graph G. Two vertices are called 
connected, if exists I > such that the vertices are l-connected. Otherwise the 
two vertices are called disconnected. By default, every vertex is considered to be 
connected to itself. 

A graph G = (V, E) is called connected, if all vertices in the graph are con- 
nected. Otherwise the graph is called disconnected. 

A graph vertex v is of degree k, deg(w) = k, iff there are k different edges 
that contain v. 

A degree partition D(G) of a graph G — (V, E) is formed by an equivalence 
relation among the graph vertices {(^1,^2) G V 2 | deg(i>i) = deg(i>2)}. 

Definition 2.5. A complement of a graph G — (V,E) is a graph G — (V,E), 
where E = {{v%, v 2 } \ {vi, V2} ^ E & vi ^ v% & V\,vi G V}. 

Theorem 2.6. // a graph is disconnected, then its complement is connected. 

Proof. Let G = {V 1 E) be a disconnected graph and v\ , v-i - two vertices of the 
graph. If v\, V2 are disconnected in the graph G, they form an edge {vi, V2} in a 
graph G and thus are connected in G. 

Suppose vi,V2 are connected in G. Since G is disconnected, exists a vertex 
vq, which is not connected to v\ in G. Hence i>o is not connected to V2 in G as 
well. Therefore the vertices vq, v\, V2 form edges vq} and {vq, V2} in G. Thus 
r>i,vo, V2 is a path on the graph G and the vertices vi, V2 are connected in G. 

Hence the graph G is connected. □ 

Definition 2.7. We shall call a graph doubly connected |L OOj if both the graph 
and its complement are connected. 

We can apply a bijection V\ V2 on the set of edges of a graph G = 
i.e., {vi,v 2 }P = {vi/3,v 2 P}, hence E x [3 = {{vif3,v 2 P} \ {v 1: v 2 } G Ei}. 

Definition 2.8. Two graphs G\ — (Vi,Ei), G2 = (V^-Eb) are called isomorphic 

if exists a bijection V\ — ► V2 such that E\tt = E2- The bijection (3 is then called 
an isomorphism. 

Thus graph isomorphism problem for two graphs G\ and G2 means finding 
whether exists an isomorphism between these two graphs. 
A permutation 7r on V is a bijection V — > V. 



Definition 2.9. A permutation it on V is called an automorphism of a graph 
G = (V,E), if En = E. 



We denote by A(G) the set of all automorphisms of a graph G. The following 
lemma is straightforward. 

Lemma 2.10. Given a graph G = (V, E) and an automorphism n £ A(G), for 
all Vi,V2 £ V i>i,t>2 are 1-connected iff Viir, V2^ are 1-connected. 

Corollary 2.11. Given a graph G — {V, E) and an automorphism n £ A(G), 
for all v E V deg(w) = deg(«7r). 

The automorphism relation a among the vertices of the graph is defined as 
follows. 

Definition 2.12. For two vertices Vx,V2 £ V of a graph G — (V,E), v\omi iff 
exists an automorphism n £ A{G) such that v\-k — V2 ■ 

Since a is an equivalence relation, it spans the set of vertices V into a set of 
equivalence classes, called automorphism partition of the vertices of a graph. 

Thus automorphism partitioning problem for a given graph is to find a par- 
tition of the graph's vertices into automorphism equivalence classes. 

We denote by P(G) the automorphism partition of a graph G. It is quite 
evident that given a graph G, P(G) = P(G). By Corollary 12. Ill if two vertices 
are in the same class of automorphism partition, they are in the same class of 
degree partition. 

Definition 2.13. Two partitions of P{V\), P{V2) of the vertices of separate 
graphs G\ — (Vi,E\) and G2 = (V2,E 2 ) are called equivalent, if exists a bi- 

jection V x V 2 such that P{V{)(3 = P(V 2 ). 

Two isomorphic graphs have equivalent automorphism partitions, the oppo- 
site is not necessarily true. 

In |RC 77] it is shown that graph automorphism partitioning problem and 
graph isomorphism problem are polynomially equivalent. As mentioned in Sec- 
tion ^ algorithm Al is a vertex classification algorithm toward automorphism 
partitioning classes. For the reasons exposed in Section Al is applied only 
on doubly connected graphs. Therefore it is necessary to state the following 
theorem. 

Theorem 2.14. The graph isomorphism problem is polynomially reducible to 
the automorphism partitioning problem of a doubly connected graph. 

Proof. Assume a polynomial automorphism partitioning algorithm is given. Con- 
sider two separate graphs G\ = (Vi, E\) and G2 — (V2, E2), \G\ \ — |G?2 1 = n > 1. 
Assume that graphs have equivalent degree partitions, otherwise the graphs are 
not isomorphic. Further assume that both G± and G2 are connected, otherwise 
consider their complements G± and G2. Select a vertex Co £ V\ with a maximal 
degree in G\, deg(co) = do, and for each vertex r £ {v £ V2 | deg(u) = do} 
form G T = {V,E T ), where V = Vi U V 2 and E T = E x U E 2 U {{ct ,t}}. G t is 



connected. The graph G T is connected. (Proof is similar to that of Theorem ^.GI l 
Therefore G r is doubly connected. 

Now for every r perform the automorphism partitioning on G T . If exists r 
such that agar, then G\ and G2 are isomorphic, otherwise the graphs are not 
isomorphic. (Lemma 12 .151 ) 

The number of performed automorphism partitionings is bounded by n, hence 
the reduction is polynomial. □ 

The following lemma certifies the correctness of the reduction algorithm de- 
scribed in the proof of Theorem 12. 141 

Lemma 2.15. The graphs G\ — {V\,E{) and G 2 — (V 2 ,E 2 ) are isomorphic iff 
exists t £ V 2 such that a^ar holds for G T . 

Proof. Suppose that exists r G V 2 such that a^ar holds for G T . Hence exists 
7r G A(G T ) such that Co 71 " = T - Vertices <7q,t are the only vertices with degree 
do + I in G T . Therefore they are the only members of an automorphism class 
in P(G T ). Hence ttt = cto- Since V\ is connected, it is possible to construct a 
sequence {cr } = lo C Ii C ... C Xn-\ = Vi, where a set X i+1 is formed 
by adding to Xi a vertex in V\ \ Xi which is 1-connected to a vertex in Xj. 
Let us prove that Vi > Xi?r C V 2 . Proof is by induction. As for induction 
basis, (To 71 " = t E E 2 . Now, for induction step, assume that Xiit C V 2 . Take 
the vertex v G Xi + \ \ Xi. Consider a vertex x G Xi, which is 1-connected to v. 
Since xit e V%, by Lemma |2.10l either vk = ctq or vir G V2. However ->(vaao), 
so vtt G V2 and Xi + i7r C V%. Now, V±ir C V2. Since n is a permutation and 
\Vi = I V2 1 , V\tt = V 2 . Furthermore, that leads to V 2 tt = V\. Since E t tt = E T , 
Eitt U E 2 tt = Ei U E 2 . Since £ 1 n£ 2 = Eiir n E 2 n = Ex n Bi7r = £ 2 n E 2 7r = 0, 
E\Tt = E2 and = -Ei- Hence 7r is an isomorphism from Gi to G2. 

Suppose Gi is isomorphic to G 2 . Then exists an isomorphism V17 = V2, 
£Ji7 = i?2- Define a permutation 7r on V\ U V 2 as 7r = 7 U 7 _1 . Now V\tt = V 2 , 
V 2 7r = V\, E\-k = E2 and i? 2 7r = Take r = cro 71 "- By construction, ttx = a®. 
Consider the graph G T . For the set of the graph's vertices, Vtt = {V\ U V^)^ = 



Vitt U V 2 tt = Vi U V 2 = V. For the set of the graph's edges, E t tt = (E 1 \JE 2 U 



{{o- ,t}})tt = E 1 Tr\JE 2 iTU{{cro,T}n} = E x I) E 2 U {{a ,r}} = E T . Hence vr is 



3 Dynamical System and Its Basic Properties 

In this section we define the dynamical system and prove basic properties of the 
system upon which the partitioning algorithm Al will depend on. 

We shall take use of some definitions and theorems from the theory of real 
analytic functions. For detailed exposition, see |KP 02j . 

Suppose a graph on m vertices is given and its adjacency matrix is H = 
(hik). The basic idea behind is to consider the graph as a physical system, 
embedded into m-dimensional space. The reason why dimension should coincide 



an automorphism of G r . Therefore o^crr. 



□ 



with the number of vertices will be explained in Section 0] Vertices of a graph 
are points that enjoy mutual attraction and repulsion forces. Every two points 
are interacted by repulsion force, whereas every two points corresponding to 1- 
connected vertices are also interacted by attraction force. The mass of any point 
is defined to be 1. 

The coordinates of m points form an to x to matrix X, where the i-th row 
of the matrix denotes the coordinates of the i-th point, Xi — (xu, Xi2, ■ ■ ■ , %i m ). 
Repulsion force Fr^ between any two points x\ and Xk is chosen to be 

/ m 

inversely proportional to the distance \\xi — Xk\\ = \ ( X U ~ x ki) 2 between the 

V 1=1 

two points, ||-FVifc|| = l/||xi — x k \\, whereas attraction force Fa,i k between the 
points connected by an edge is chosen to be directly proportional, ||.Fajfc|| = 
\\xi - x k \\. 

In principle, it is possible to define the forces in a different way, say, H-Fr-jfeH = 
l/||a;i — Xfc|| 2 or ||-FYifc|| = 1. In this paper, we did not attempt to investigate 
the impact of choosing particular attraction or repulsion force, however it seems 
that attraction and repulsion forces should be linearly independent. 

So attraction force acting on a point Xi due to its interaction with a point 
x k is 

Fa ik = -h ik (xi - x k ), (1) 

whereas repulsion force 

PW = - 

\\Xi - x k \\ 2 



Fr lk = .. Xi X \ 2 . (2) 



The total force acting on Xi from x k is 

1 

'\Xi ~ XjcW 

Hence Fi k = —F k i and Fa — 0. So the total force acting on Xi is 



F ik = Fr ik + Fa ik = (x t - x k ) ( — - h ik ). (3) 



mm f \ \ 

* = E *•* = E - »*) (r^p - h *) ■ ( 4 ) 

In our setting, 

d 2 x t {t) 
dt 2 

hence the dynamics is described by the following autonomous system of differ- 
ential equations; 

i,i j 4 = g ( Xij - x kj ) ( ^^p - ^) • (5) 



Let M. be a set of to x to real matrices, whereas M + C M. a set of matrices 
where no two rows in the matrix are the same. The right hand side of J5J may be 



regarded as a function Fh ■ Ai + — > M. Hence the system J5J may be rewritten 

as 

X" = F H (X). (6) 
By making a standard reduction to first-order differential equations, we obtain 

X' = Y 

y = f h (x). (7) 

Given B G Ai + , V G Ai, initial conditions for JJJ) are specified as 

(8) 

The system (JJJ is in fact a Hamiltonian system, that is, 

dU 



* (to) = B 
Y(t ) = V. 



x ij 



y'ij 

i 



dyij 
dU 



(9) 



Oxij 



where yij are elements of Y and 

H = E k (Y) + E p (X) , (fO) 



1 m 

my) y 2 r* . ( n ) 

r,s — 1 

^ 711 

E P (X) = -- (log ||x r - a; s || 2 -/i rs H^-^ || 2 ). (12) 

r,s— 1 
r<s 

Functions -Efe(F) and E p (X) are frequently referred as kinetic and potential 
energy, respectively. The system {JJ being an autonomous Hamiltonian system, 
energy conservation law f |HS 74| . p. 292; |TP 63) . p. 475) is applicable, yielding 

E k (Y) + E p (X) = E , (13) 

where Eq is a constant value; using initial conditions JHJl one can compute Eq = 
E k (V)+E p (B). 

We would like to prove that H7I8J1 is a dynamical system, that is, coordinates 
of points are defined for any time moment igl. 

Let Z H (t, B, V) = (X H (t, B, V), Y H (t, B, V)) be a solution of 03). The right 
hand side of @ may be regarded as a function Gh{X,Y) — (Y, Fh(X)), Gr ■ 
Ai + x Ai — > Ai x M.. Each scalar component of Gh is in essence a multivariate 
rational function, hence all of them are real analytic functions onM + x M.. So 
each component of Gh is continuously differentiable on Ai + x M C |KP 02j . p. 30). 
Therefore G H is locally Lipschitz on M + f |HS 74| . p.163). So for any B £ M + 
V" G X exists e > such that C |HS 74j . pp.162-163) 



1) Zn{t, B, V) exists on some open interval (fi, £2), where t\ = to — e, = to+e; 

2) Z H (t,B,V) is unique. 

The local solution defined on (tijfa) can be extended to a unique solution 
defined on maximal open interval (a, (3), where a < t\, /3 > ia ( |HS 74j . p. 171). 
That also implies that X H (t, B, V) £ M + for all t £ (a, 0). 

Lemma 3.1. Suppose that Zn(t, B,V) is defined on interval [to, [3) which can- 
not be extended to the right. Suppose that (3 < +00. Then 

i) Zn(t, B,V) is bounded on interval [to, (3), that is, 3D > \ft, to < t < (3, 
\\Z H (t,B,V)\\<D; 

ii) Xn(t, B,V) stays sufficiently far from the boundary of A4 + , that is, exists 
L>0Vt, t Q <t</3, V(r, s), r^s, \\x r (t) - x s {t)\\ > L. 

Proof. The proof strongly relies on the constraints imposed by energy conser- 
vation law (| 1 31) - The proof consists of two parts. In the first part we prove 
proposition i) and in the second part - proposition ii). 

Note that if r ^ s \\x r (t) - x s (t)\\ > for any t £ [to, (3). 

/ m 

1) Let N(t) = \\X H (t,B,V)\\ = J E Xrs(t) 2 - N(t) is defined on [t ,/3). 



Due to Cauchy inequality ( B B 71) . p. 2) 

\N'(t)\ = 



E x rs (t)x' rs (t) 

r,s— 1 



N(t) 



<* E K s {t)) 2 = V^iD- 



r,s—l 



On the other hand, due to general means inequality ( BB 71 , pp. 16-17), 

m 

E k {t)=Eo-E p (t)=Eo + \ E (log\\x r (t) - x s (t)\\ 2 - h rs \\x r (t) - x s (t)\\ 2 ) < 

r,s=l 
r<s 

m m 

E + E log\\x r (t)-x s (t)\\=E + \og II \\x r (t) - x s {t)\\ < E + 

r,s— 1 r,s— 1 

r<s r<s 



log T, \\x r (t)-x s (t)\\ 



r,s=l 



3 







te s(m^r£i K(t)-*.(t)||J <E 

! ^iog (s^zij r E i (iK(t)ll + IM*)||)) = 



m(m— 1) 
2 

m(m— 1) 
2 



log ^EIM*)II ^ 



m(m—l) i 
2 bg 



E 



'E IM*)|| 2 = E 

s=l 



log (^N(t)\. Let K = E 



m(m — 1) 



log-?-. So 



< E k (t) < K + log N(t) <K+\ + logiV(t). (14) 



We have obtained a differential inequality 



\N'(t)\ < y/2K + 1 + m(m - 1) log AT(t). (15) 



The corresponding initial condition is N(to) = \\B\\, where = , / E {B rs ) 2 . 

y r,s=l 

Hence 

N'(t) < V / 2K + 1 + m{m - l)\ogN(t), N(t ) = \\B\\. (16) 

Consider initial value problem 



N[(t) = ^2K + 1 + m(m - 1) log JV x (t), A^o) - \\B\\. (17) 

Due to ijTijl. 2AT + 1 + m(m — 1) log ||B|| > 0, so (fTTf satisfies Lipschitz condition 
and has a unique solution on (to — e,t + e), e > 0. The solution of (|T7|l is 
characterized by equation 

t = t + I dZ (18) 

J y/2K + 1 + m(m - l)logz 

| | 

As ATi — > +oo, the definite integral diverges and t — > +oo. Therefore it is possible 
to extend the solution of l|17fl to interval (to — e, +oo). This implies that Ni(t) 
is bounded on interval [to, (3). 

Due to |LL 69j . pp.7-8, A^t) > N(t) on [t , /?). Note that N(t) > 0. So N(t) 
is bounded on [to, (3). But then Xu(t) is bounded on [to,/3) a s well. 

Boundedness of Xjj(t) implies that exists P > such that —E p (t) < P. So 
< E k (t) = Eo- E. p (t) < E + P. Therefore Y H (t) = v / 2E k (t) is bounded. 
Since both -Xjy(t) and Yu(t) are bounded, Zn(t) is bounded. 

2) Let d(t) = min ||2v(t) — a; s (t)||. Since Xjj(t) is bounded, exists Z > such 

that for all r,s and £ G [t ,/3) \\x r (t) - x s (t)\\ < Z. Set L = ^- m( l _1> +i e --Eo_ 
Due to energy equation (| 1 3|) - 

m 

= § E (log \\Xr(t) ~X s {t)\\ 2 -h rs \\ Xr (t) -X s (t)\\ 2 ) < 



r,s=l 



E loglkrW - x s (t)\\ < log d(t) + (2fc2r!l _ i)io g Z. So logd(t) > -S - 

r.s — 1 

r<s 

^mCm-i) _ ^jog^ and > L So || Xr (i)_a; a (f)|| > L. □ 

Lemma 3.2. Suppose that Zn(t, B,V) is defined on interval (a, to] which can- 
not be extended to the left. Suppose that a > — oo. Then 

i) Zn(t, B,V) is bounded on interval (a, to], that is, 3D > Vt, a < t < to, 
\\X H (t,B,V)\\<D; 

ii) Xjj (t, B, V) stays sufficiently far from the boundary of Ai + , that is, exists 
L > Vt, a < t < t , V(r,s), r s, \\x r (t) - x s (t)\\ > L. 



Proof. Proof is essentially the same as above. We obtain inequality l|15p. Now 
N'(t) > -y/2K +l + m(m- l)logiV(i), N(t ) = \\B\\. Through substitution 
s = —t, sq = —to, we define M(s) = N(—s) = N(t), where M(s) is defined on 
interval [s , -a ). So Af'(s) = -N'(t). We ha ve 

-M'(s) > -^2K + 1 + m(m - 1) log M(s), M(s ) = ||B||, yielding M'(s) < 
y/2K+l + m(m — l)logJVf(s), M(so) = \\B\\. We have actually obtained the 
inequality (|16(l . Therefore M(s) is bounded on [sq, —a) and iV(t) is bounded on 
(a, to]. So Xu(t) is bounded on (a, to]. The rest of the proof is the same as for 
the corresponding parts of Lemma l3.ll □ 

The previous two lemmas show that no collisions are possible and no point 
may reach infinity in a finite time. This is different from the classical Newtonian 
n-body problem, where both types of singularities are possible |X 92| . 

Theorem 3.3. The solution Zn(t, B ,V) is defined on interval (— oo, +oo). 

Proof. Let us assume from the contrary that the maximal interval of the function 
Z H {t) = (X H (t),Y H {t)) is (a,/3), where j3 < +oc. 

Since @ is autonomous, as t — » +/3 either Zn(t) tends to the boundary 
of M+ x M or \\Z H (t)\\ -> +oo, or both f |HS 74 |. pp.171-172). However, due 
to Lemma 13.11 any of the two options is not possible, which is a contradiction. 
Hence j3 = +oo. 

For similar reason, due to Lemma 13.21 a — — oo. □ 

As remarked below the equation (|13|l . Fh{X) is real analytic on A4 + . So a 
consequence of Cauchy-Kowalewskaya theorem (special case for ordinary differ- 
ential equations; |KP 02j . p.42; |TP 63j . pp. 555-556) is that for any B E M+, 
V e M, the solution of'CEJ Z H (t,B,V) = (X H {t,B,V),X' H (t,B,V)) is real 
analytic at to. Furthermore, the consequence of Theorem 13.31 is that Vf G R 
Xff{t, B, V) G A4 + . Therefore, initial conditions may be stated at any point 
(ii, Zjj(t\, B,V)), so Cauchy-Kowalewskaya theorem may be reapplied to state 
that Zn(t, B, V) is real analytic on R. 

So Xn(t, B,V) is real analytic on R. Hence for each t\ G M x ij{t) ma y 

oo 

be expanded into convergent power series {ayn(^i)} = ^2 a-ijn(t — ti) n , where 

n=0 

o-ijn — — — j — ■ The series converge on some interval (t\ — e, t\ + e), e > 0. 

The power series can be compared by comparing their coefficients. Now the 
possibility to compare the power series aij n instead of the functions Xij is asserted 
by the following theorem. 

Theorem 3.4. The power series {a,ij n (ti)} and {a rsn (ti)} are equal iff the func- 
tions Xij(t) and x rs (t) are equal on R. 

Proof. If the functions x%j(t) and x rs (t) are equal, then their expansion into 
power series is the same. 

Suppose the power series {aij n (ti)} and {a rsn (ti)} are equal. Both power 
series are convergent on some interval (t\ — e,t\ + e). So the functions X{j(t) 



and x rs (t) are equal on (t\ — s,t% + e). However as stated in |KP 02) . p. 14, if 
two functions are real analytic on an open interval U and are equal on an open 
interval V C U then these functions are equal on U. Therefore Xij(t) and x rs (t) 
are equal on R. □ 

Let us now consider some of the combinatorial properties of the differential 
equation ©. 

Let fn{B,t) = Xn(t, B,Vo). So /#(£?, t) is a real matrix whose rows are 
coordinates of m points of the physical system at the time moment i G 1, with 
initial configuration of points being B. Let P be an arbitrary mxm permutation 
matrix. Relabeling of the coordinates of the points does not have impact on the 
evolution of the system, therefore it is quite straightforward that 



f H (BP,t) = f H (B,t)P. 



(19) 



Quite similarly, relabeling of the points (and the vertices of the graph H , respec- 
tively) does not have effect on the evolution of our physical system, therefore 



f H (PB,t)=Pf P T HP (B,t). 



(20) 



Equations (|19I20() play key role in the partitioning algorithm introduced in 
the next section. 



4 Partitioning Algorithm 

Due to Theorem 12. 141 graph isomorphism problem is reducible to the automor- 
phism partitioning of doubly connected graphs. In this section, a partitioning 
algorithm Al is described which is applied for this type of graphs. 

Let us consider the system (JHJ), which defines the dynamics of m points 
corresponding to a graph vertices. It is aimed that future coordinates of the 
points are characterized only by dynamics of the system and are not dependent 
from their initial configuration. The system is embedded into m-dimensional 
space exactly for this purpose; initial coordinates of the points are characterized 
by an m x m identity matrix I mi that is, the coordinates vectors of the points 
initially are 



Xi 


= (1,0,0,. 


.,0,0) 


X2 


= (o,i,o,. 


.,0,0) 


X3 


= (0,0,1,. 


.,0,0) 


*£m- 


-i = (0,0,0,. 


.,1,0) 




= (0,0,0,. 


.,0,1) 



The initial velocity of the points is 0. 




Hence we consider a case where the coordinates of the points initially form 
the identity matrix / and the initial velocity of the points is 0. That corresponds 
to a system 

X" = F H (X) 

X(0) = I (21) 
X'(0)=O, 

where O is zero matrix. 

Let A be a permutation matrix corresponding to an automorphism of the 
graph H. Hence A T HA = H . Recall equations (|19l 120ft which express some of 
the properties of the function fn, characterizing the future coordinates of the 
system depending on the initial configuration. By l |T5jl and i(2Ul) . /#(!, to) = 
f H (A T IA,t ) = A T f AHAT (I,t )A = A T f H {I,t Q )A. That is, 



(22) 



We say that two rows r\ and r-i of a matrix are permutation equivalent, if exists a 
permutation matrix P, such that t\P = r 2 . Let Xji(t) = Xn(t, I, O) = /h(/, t) 
be the solution of 1)21(1 . A straightforward corollary of the equation \\2'2\ is the 
following theorem. 

Theorem 4.1. If two vertices of a graph H are automorphism equivalent, then 
the corresponding two rows (and columns) in Xu{t) are permutation equivalent. 

The last theorem reveals the idea of the vertices partitioning algorithm - com- 
parison of those m? functions which are the solutions of the system (|21l) . 
Two vertices are in the same class of partitioning iff the corresponding rows of 
Xn(t) are permutation equivalent. 

One method to compare the functions would be numerical computation of 
Xn(to) at some time moment t$. Such an algorithm would be relatively fast. 
However, two different functions may as well be of the same value at to- It is also 
not quite clear which is the precision required for such computation, to separate 
two values, say, x\^{to) and £2,3(to), which could be very close to each other 
and yet not equal. 

However, due to Theorem 13.41 instead of comparing the functions Xij and 
x rs directly it is possible to compute Xn(t) as formal power series and then 
compare coefficients of the first terms of the series. 

°° xif (0) 

Hence \/(i,j), 1 < i,j < m, Xij(t) — ^2 a>ijnt n , where Oij n — — — -j — ■ Let 

€L\Xn • • • ftlmi 

A(n)= \ \,soX H (t)= £ A(n)t» 

^ml)i ■ ■ ■ drnmn 

Let us solve the problem l|21|l in terms of formal power series. We denote 
a ij(2n) as a^j, A(2n) as A n , and seek for the solution in the form Xij(t) — dij = 



E Q-ijt 271 1 thinking of Xij as series with unknown coefficients. Due to the initial 

n=0 

conditions, Aq = I. 

Consider the equation 10). Let us compute the s-th term of the left and the 
right side of J5J. As for the s-th term of the series x"j, 



[4] s = 2( S + l)(2 S + l)a^ 1 . 
Recall that for two power series 



n=0 



and if ap 0, b s 



E a nt U 



n=0 



ra=0 
1" 



(23) 
(24) 



d=0 



may be computed by recurrence 



1 s 1 

b s = a c b s ^ Cl where b = — . 

a ' a Q 

c=l 



(25) 



Now 



- ^f] c - E E(^ - - 47 d )- (26) 

2 = 1 d=0 

So [|| acs - a? fc || 2 ] = 2(1 - S ik ) and r 8 ik = [\\xi - a; fc || _2 ] s , i ^ fc, is computed by 
recurrence 



^Err[ll^-^H 2 ] C , wherer° fc = l (27) 

c=l 

We define = 0. Therefore using the right side of ©, we obtain 

/ m s \ m 

[41- EE^-W -EM4-«U ( 2 8) 



\k=l p=0 



k=l 



In complete analogy to A n , we introduce matrices R n — (r™ ) and define 1 as 
an to x m matrix where Vi, j ljj = 1. So equations (|23|l and (|28|l give a doubly 
recurrent formula to compute the s-th term of the series x^, 



n s+1 - 1 

U ij 2(s+l)(2s+l) 



E J2(a P lj -a p k3 )r s ik p 



k=l p=0 



E ft ifc( a 



k=l 



(29) 



c=l 1=1 d=0 



where Ao ~ I and i?o = |(1 — 

Our solution implies that the function's iy (<) derivatives of odd order at 
t = vanish, x^ q+1] (0) = 0. The formula igSJ is in essence an algorithm to 
compute A n . 

Algorithm Al. Now the vertex partitioning algorithm Al is as follows: 



1) compute A s for < s < m 2 ; 

2) place two vertices i and j into one partition class if and only if Vs, < s < m 2 , 
the i-th and the j-th row of A s are permutation equivalent. 

The arguments why checking the first m 2 elements of A s should be sufficient 
to distinguish between two functions that are not equal are discussed further in 
the section. 

Now we are going to prove that Al is a polynomial-time algorithm. For that 
purpose, we state the following lemma. 

Lemma 4.2. If a, 6 are inteqers such that a > 0, b > then - — ^— — ^— is an 

~ a\b\b\ 

even integer. 

Proof. The number - — rrr-rr- is a multinomial coefficient and therefore is an 
a! 6! 6! 
(26)! 2(26 - 1)! 

integer. If a = 0, , , , , = -77- -— , so it is even. Let q be a positive integer and 

6! 6! 6! (6 — I)! 

(g + 26)! . _. (g + 1 + 26)! (g + 26)! 2(g + 26)! 

suppose — j— — is even, bmcc — 



9 !6!6! ' (q + l)\b\b\ qlblbl (q + 1)! 6! (6 — I)! ' 

(g + 1 + 26)! . „ , . . . (a + 26)! . 

— ,,,,,, is even. Hence by induction — . , . , . is even. □ 

(g + l)!6!6! a! 6! 6! 

Theorem 4.3. The algorithm Al is polynomial in time with respect to m. 

Proof. It is straightforward to check that the number of additions and multi- 
plications necessary to perform to compute A s is polynomial with respect to s. 
Since s < m 2 , the number of arithmetic operations is polynomial with respect 
to m. 

Therefore to show that Al is polynomial it is sufficient to prove that each 
arithmetical operation involving af • and rfj can be performed in polynomial 
time with respect to s. The numbers and rf^ are rational. We may assume 
that for each s a?- and rf„- are calculated to their reduced fraction form. If 
a rational number is negative, we assume that its numerator is negative. Let 

af, = -r— and rf, = -r— be those reduced fractions. Hence we have to show 
d is } d 2s 

that the logarithms of |ni s |, \n2s\, d± s , d2 S , i-e., the lengths of numerators and 
denominators, are bounded from above by some polynomial P(s). 

Let c4 J = 2 s (2s)! o?- and pf k = 2 s+1 (2s)! rf k . So a% and p° ik are integers. The 
equations l|29|) imply that 



a s+1 



EE Jk-W ) -2X)m«?j-«u (30) 

fc=lp=0 V "/ / fe = l 

,+i _ 1 y ,-c+i y y W8 + 1))! (a| - a d kl )(atT d - ^j d ) 

2 hi (2rf)!(2(c-d))!(2( s -c+l))! ' 1 > 

Assume that Vt, < t < s , a^- and p\ k are integers. So by (J20J), is an 

(2s)! 

integer. Consider the equation 1)31(1 . The numbers . . -. —. . . -. —. 777 

(2d). (2{c — a)). (2(s — cj). 



are in fact multinomial coefficients and therefore are integers. If c is an odd 
integer, c = 2q + 1, 

(2d)\(2(c-d))\(2(s-c+l))\ 
^ 2 ^ (2( S + l)y.(a^ d -ali d )(af l -a d kl ) 



^ (2d)!(2(c-d))!(2(«-c+l))! 



If c is even, c = 2g, 



^ (2( g + i))!(^-^)K r rf - a - rf ) 

(2d)!(2(c-d))!(2( S -c+l))! 



_J^ {2{ S + l))\{a d l ~a d l ){a^ d -al- d ) (2(s + 1))! (a* ~ ^) 2 



d=0 



(2d)!(2(c-d))!(2(s-c + l))! c!c! (2(s- c+ 1))! 



(2(s + I))! 

By Lemma f4. 21 , , , , — is even, therefore in both cases 

— c! c! (2(s — c + 1))! 

(2( g + l))!(^-a^)(ar d -^7 d ) fa even go , » +1 
2- (2d)! (2(c - d))l (2( s - c + 1))! 18 6Ven - S ° by ^ 



is even. So by (fTTTf) p s ik is an integer. 



However, a tj = -fi- and = -^-y. So d ls < 2'(2s)l, d 2s < 



d=0 

Hence by induction Vs > ah and pf k are integers. 

anr | r s . Pik 

2 s {2s)\ lk 2 s + 1 {2s)\ 

2 s+1 (2s)! and logdi s < ci (slogs), logd 2s < c 2 (slogs). 

On the other hand, afj is generated by real analytic function Xij (t) , whereas 

r s ik is generated by a function r^. (t) = - — — , which is real analytic 

for the same reason as Xij(t). Therefore ( jKP 02j . p. 15) exist constants C > 

C C C2 s (2s)' 
and R > such that af„ | < — and \rfA < -—. Hence afJ < and 

C 2 s+1 f2s)! 

< ^ '-■ So log|ni s | < c 3 (slogs) and \og\n 2s \ < c 4 (slogs). 

Theorem is proved. □ 

By Theorem l4.ll if two vertices are automorphism equivalent, the algorithm 
Al places them into one partition class. It is however an open question what are 
the types of graphs the algorithm Al produces an automorphism partition for. 

The condition that a graph must be doubly connected is necessary. Suppose 
a graph G consists of two non-isomorphic strongly regular graphs of the same 
parameters. The graph G is not connected, hence it is not doubly connected. 
As experiments show, in this case, algorithm Al does not give automorphism 
partition of G (and G as well) , which happens for essentially the same reason as 
why Gudkov-Nussinov algorithm fails |S.IC 03j . 

Currently we do not have any additional satisfactory explanation, yet the 
following intuitive argument might be helpful. If a graph is not connected it con- 
sists of several connected subgraphs. So intuitively, there is no enough interaction 



between these subgraphs for algorithm Al to work in this case. The situation for 
complement of a disconnected graph is dually the same, now both repulsion and 
attraction forces are in effect between every two vertices each belonging to one 
of the corresponding subgraphs. Again, the interaction among the subgraphs is 
too homogeneous. 

The algorithm Al only computes the first to 2 coefficients of the series a,ij. 
It is not formally proved whether this is sufficient. However, below a symbolic 
algorithm Al' devised from Al is presented, where first m 2 are proved to be 
sufficient to distinguish between two different functions. 

s 

Let afj be the s-th partial sum of ay, i.e., af- = Yl a ijt 2s - So a ^ = a°-. 

71=0 

Let Z s be the set of s-th partial sums of formal power series and define the 
convolution and deconvolution operations *,0:Z s xZ s — > Z s , where ay * by, 
called convolution, is the s-th partial sum of the multiplication ciijbij and ajj06jj, 
called deconvolution, is the s-th partial sum of the polynomial division of ay by 
bij. 

It is straightforward to verify that (J5J and 129|) imply that 

t t 

«;/' 4 + //(E (K- - O s tk - K - (32) 
° ° %\ 

m 

where S? k = ||af - a s k \\ 2 = £ ((aj, - a s kl ) * (a s tl - a s kl )). 

i=i 

Now if ai r 7^ cijs then 3k, such that af r ^ a^ s . To evaluate k, let us trans- 
form H^2(l into the following symbolical algorithm. Let U, P be letter sequences 
{ao; &i, . . .} and {/3o, 01, ■ ■ ■}, respectively Consider matrices C s and D s with 

elements in U and P, respectively. Let = < a °' 1 , . . For any s, we compute 

D s in the following way. We define ordering on N 2 ; (k, I) < (i, j) iff (k = i & I < j) 
or k < i. 

1. FORALL let Py = (cfc.cfe), 

2. FORALL i, j, sort Py as pairs; 

3. n = 0; 

4. FOR i = l TO rn 

FOR j = 1 TO m 

IF 3fc, Z such that (k, I) < and Pu 

u ij ~ u kl 

ELSE 

d^j — fin 
n = n + 1; 

Now C" 5 "*" 1 is computed as follows. 

1. FORALL let 

T y = c?i> (c?j, d? 2 , /i <2 )> ■ • • > ( c ?j' c mj> a vector 

of 4-tuples; 



• • • > (C c fcm)] a vector of P airs ; 



= P i7 - THEN 



2. FORALL sort as tuples; 

3. n = 0; 

4. FOR i = 1 TO m 

FOR i = 1 TO m 

IF 3k, I such that (k, I) < and T ki = T l3 THEN 

„s _ s 
L ij — L kl 

ELSE 

c ij = a n 

n = n + 1; 

Similarly as in Theorem l4.ll if two vertices i, j are automorphism equivalent 
then the corresponding two rows (and columns) in C s are permutation equiva- 
lent. 

Algorithm Al\ 

1) compute C s until C s = C s+1 ; 

2) place two vertices i and j into one partition class if and only if Vs, < s < m 2 , 
the i-th and the j-th row of C s are permutation equivalent. 

The algorithm Al' is polynomial and now checking the first m 2 elements is 
clearly sufficient: 

Theorem 4.4. The algorithm Al' is polynomial and the number of steps is 
bounded by m 2 . If X{ r (t) and Xj S (t) are not equal, then 3k, < k < m 2 , such 
that c\ r ^ c k JS . 

Proof. The symbolic algorithm is actually terminated as soon as for some s 
C s = C s+1 . Indeed, in that case for any t > s C* = C s . Further, if for some s 
c tj 7^ c ki then c^ +1 7^ c^l . It can be further noted that the algorithm performs 
no more than m 2 steps, since there can be no more than m 2 mutually different 
elements in C s and in each step, letters from £ are assigned in the same order. 

In a way, the symbolic algorithm mimics the arithmetical expressions of (|32|) . 
So if Xi r (t) ^ Xj S (t) then also afj ^ a| ; , which in turn implies ^ c s kl . But 
then exists t < m 2 such that c\j ^ c^. □ 

The algorithm Al' essentially covers the algorithm Al. Indeed, if for some 
graph Al gives an automorphism partition, so does Al'. Hence the discussed 
dynamical algorithm is covered by a symbolical (combinatorial) algorithm, which 
in essence means that the chaotic behavior of the point-masses in the presented 
dynamical system does not help too much. 

As experiments show, contrary to the approach presented in |GN 02) . the 
algorithms Al and Al' may be used to distinguish between two strongly regular 
graphs with the same parameters. (By Theorem 12 .141 graph isomorphism prob- 
lem is reduced to the automorphism partitioning problem of a doubly connected 
graph.) However, exists a counterexample where Al' (and therefore also Al) fail. 

Example 4-5- Apply the algorithm Al' to the pair of graphs A25 and B25 in 
|M 78| . which are not isomorphic. Automorphism partition of the resulting dou- 
bly connected graph is not achieved and hence the solution to the graph isomor- 
phism problem is not achieved. 



5 Yes-No-Don't Know Algorithm 

Algorithms Al and Al' may give two answers to the question whether two graphs 
are isomorphic, i.e, No or Don't Know. If the answer is No, then with all certainty 
the graphs are not isomorphic. 

However, it is possible to give an algorithm A2, which gives three answers: 
Yes (in that case the new algorithm supplies isomorphism 7 between the two 
graphs), No (with all certainty the graphs are not isomorphic) or Don't Know. 

First, let us describe an algorithm which reduces the computation of 7 to the 
automorphism partitioning of a doubly connected graph. A degree of a vertex a 
of a graph G will be denoted as deg(<r, G). 

Verification Algorithm. Suppose two separate graphs G\ = (Vi,Ei) and 
Gi = (Vi^E-i) are given, \G\\ = \Gi\ = n > \. Without affecting generality it 
is possible to assume that both G\ and G2 are connected and have equivalent 
degree partitions. Apply the reduction algorithm of Theorem 12 .141 to construct 



Now apply automorphism partitioning algorithm on G r , which is doubly 
connected. 

If no r exists, such that ao and r are in the same automorphism class of G T , 
then by Lemma 12.151 the graphs are not isomorphic. 

Now suppose exists r, such that ao and r are automorphism equivalent in 
G T . Then by Lemma 12.151 the graphs are isomorphic. Let tq = t, Zq = G T , 
whereas the set of edges of Z a denoted as i?o = E\ UE 2 U {{ctq, to}} - The vertex 
Co has maximal degree in Gi, and so does To, deg(<7o,Gi) = deg(ro,G2) = do- 
Hence ao and To are the only vertices with degree do + 1 in Zq and therefore 
they are the only vertices in their class of P(Zq). The process which leads to the 
computation of explicit isomorphism from G\ to G2 is as follows. 

Step 1. Let o\ € V\ be a vertex with the largest degree in Gi, excluding 
(Jo. For every automorphism 7r 6 A(Zo), such that uqt: — To, V\ir — V%. Hence 
exists T\ G Vq. such that a\ and t\ are in the same class of P{Zq). So exists an 
automorphism m G A(Zo) such that Vi, < i < 1, a^i = n and t^i — &i. 
Construct Z\ = (V\ U V2, Ri) from Z by setting Ri = R U {{cto, n}, {cti, r }}. 
So deg(cr , Zi) = deg(r ,Zi) > deg(cri, Zi) = deg(n,Zi). Since R tt 1 = R , 
tt\ £ A(Z\). Apply automorphism partitioning algorithm on Z\. 

Step 2. Now let (T2 G V\ be a vertex with the largest degree in G±, excluding 
ao,cri. Again, since V\-k\ — V2, exists T2 £ V 2 such that 02 and T2 are in the 
same class of P(Z%). Hence exists an automorphism TT2 G A(Zi) such that Vi, 
< i < 2, 0^2 = Ti and Tj7T2 = cr^. Construct Z 2 = (Vi U V2,i?2) from Z\ 
by setting R 2 = i?i U {{cr^ra-J < i < 2}. So deg(a 0) Z 2 ) = deg(r ,Z 2 ) > 
deg(cri,Z 2 ) = deg(ri,Z 2 ) > deg(a 2 ,Z 2 ) = deg(r 2 ,Z 2 ). Since i?i7r 2 = Ri, tt 2 € 
A(Zi2). Apply automorphism partitioning algorithm on Z2. 

Before the j-th step, we have obtained a graph Zj-i such that Vi, < i < 
the vertex Oi 6 V\ is 1-connected to vertices {t^ s V2 < k < j — 1 — i}. 
Given deg(f7i,Gi) = deg(r 4 ,G 2 ) = d i: for all i, < i < j - 1, deg(a;, ^-1) = 
deg(ri,2j_i) = ^ + j - i. Since Vi, 1 < z < j — 1, <ii-i > d 4 , deg(cr 4 _i, = 




deg(rj_i, Zj-i) > deg(ai, Zj_\) = deg(rj, Zj-i). Furthermore, Vi Oi and n are 
the only vertices with a given degree in so {<7j,Ti} £ P(Zj_i) and exists 

automorphism 7Tj_i G A(Zj_i) such that Vi, < i < j — 1, <7,7rj_i = Tj and 
ViTTj^i = Vi- Let a, G Vi be a vertex with the largest degree in G\, excluding 
{oi < i < j — 1}. Since V\Wj-i = V2, exists tj G V2 such that aj and t,- are in 
the same class of P(Zj_\). So exists an automorphism ttj G A(Zj-\) such that 
Vi, < i < j, diTTj = Ti and Tiirj — <Ti. Construct Zj = {V\ U V2, -Rj) from 
by setting i?j = U {{(Ti, Tj^i} < i < j}. SoVi, 1 < i < j, deg(ai_i, Zj) = 
deg(ri_i,Z/) > deg(cr i ,Z :) ) = degfa, Z,-). Since Rj-iitj = Rj-i, nj G 

The process is continued until there are no two different vertices in V\ which 
are in the same class of P(Z p ), where the number of steps p < n — 2. This bound 
is achieved if a complete graph K n is tested for isomorphism to a copy of itself. 
The computed isomorphism 7 is identified by P{Z p ) because every class of Z p 
has exactly two vertices, one in V\ and other in V2- 

Algorithm A2. For every j, < j < n — 2, Zj is doubly connected. To prove 
that, it is sufficient to consider Z„_ 2 constructed from G\ = G2 = K n . In 
Z n -2, 0"n-i S V\ is 1-connected to every vertex in V2, whereas t„_i G V2 is 1- 
connected to every vertex in V\. Hence Z„_ 2 is connected. Therefore verification 
algorithm above can be transformed to the Yes-No-Don't Know algorithm A2 
by using partitioning algorithm Al instead of generic automorphism partitioning 
algorithm. If after some step j the process could not be continued for any reason 
(for example, cjj and Tj are not in the same class of partition attained by applying 
Al to Zj), the answer is Don't Know. After the last step, the computed bijection 
7 is tested for being an isomorphism. If 7 is an isomorphism, the answer is surely 
Yes, otherwise it is Don't Know. 
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